Theory of oxygen K-edge x-ray absorption spectra of cuprates 
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The dynamical mean-field theory of the three-band model of copper-oxide superconductors is used 
to calculate the doping dependence of the intensity of the oxygen _ft'-edge x-ray absorption spectra 
of high- Tc copper-oxide superconductors. The model is found not to reproduce the results of a 
recent experiment, suggesting that at sufficiently high doping the physics beyond the conventional 
three-band model becomes important. 
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I. INTRODUCTION 

The important structural unit of the high temperature 
copper-oxide superconductors is the "Cu02 plane" and 
the important electronic orbitals are believed to be the 
Cu 3(ij.2_j^2 and O 2po- states. High-Tc superconductors 
are formed by doping (adding carriers to) an insulating 
parent compound. The nominal electronic configuration 
of the insulating parent compound is Cu 3d® (one hole 
in d.j.2_y2 band) O 2p^ (oxygen p-band completely full), 
although of course Cu-0 hybridization means that some 
fraction of the hole delocalizes onto the O. However, the 
strong correlations characteristic of the Cu d orbitals im- 
plies that the energy required to add a second hole on the 
Cu site is very large, U ~ 9eV— so that at low doping 
levels, doped holes go primarily onto the oxygen site, but 
are strongly coupled via an exchange interaction to the 
holes on the Cu sites, forming "Zhang-Rice singlets".'^ 
While this picture is qualitative, being strictly valid only 
in the strong coupling limit and low doping limits, it has 
provided a useful guide for thinking about the materi- 
als. An important open question concerns the doping at 
which this picture breaks down. 

Oxygen i^-edge x-ray absorption spectra provides an 
interesting test of the Zhang- Rice picture, and more gen- 
erally of our understanding of the three-band model. In 
these experiments absorption of an incident photon pro- 
motes an electron from the oxygen Is shell to an unoc- 
cupied oxygen orbital; if the photon energy is appropri- 
ately tuned, the final states are in the energy range of 
the O 2p manifold and the spectrum reveals the energy 
distribution of the unoccupied O 2p states (modified by 
the excitonic core-hole/excited state interactions). Ex- 
periments performed in the early 1990s^ on insulating 
and lightly-doped samples La2_a;Sra;Cu04±5 (LSCO) re- 
vealed a two-peak structure, with a higher energy peak 
visible in both insulating and hole-doped materials and a 
lower energy peak, visible only in hole doped compounds. 
The integrated spectral weight of the lower energy fea- 
ture was found to increase linearly with doping (at low 
hole-dopings) whereas the integrated spectral weight of 
the higher energy feature was found to decrease4i^ For 
this reason the higher energy peak was interpreted as the 
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FIG. 1: Relative intensity of the peak associated with the 
Zhang-Rice band of O K-eA^e x-ray absorption spectra v.s. 
doping. Theplus-sign shows experiments on LSCO by Chen 
et al (Ref. [j) and the "x" -shaped sign shows experiments 
on LSCO by Peets et al (Ref [a). The diamonds, circles, 
squares and triangles are theoretically calculated integrated 
spectral weight of unoccupied oxygen states, normalized such 
that their values at 0.22 hole doping are the same. Dashed 
lines are only guides to eyes. Parameters: U — 9eV, tpd = 
l.GeV, T = O.leV^, tpp and A = — as shown in the 
legends. 



"upper Hubbard band" , while the lower energy feature 
was interpreted as the "Zhang-Rice band" . 

A recent experimenli^ has extended the measurements 
to a wider range of compounds and in particular to a 
wider range of hole dopings. Fig. [1] summarizes the data, 
showing the recent measurements as "x" -shaped crosses 
with errorbars and the previous data'* as plus signs. A re- 
markable finding of the recent measurements is that the 
oxygen weight in the feature interpreted as the Zhang- 
Rice band saturates as doping is increased beyond the op- 
timal doping value which maximizes the superconducting 
transition temperature Tc- Ref . interpreted this finding 
as a breakdown of the Zhang-Rice singlet picture. 

Motivated by these experiments, in this paper we use 
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a more modern theoretical technique, namely single-site 
dynamical mean-field theoryli^ to compute the oxygen 
iiT-edge x-ray absorption spectra implied by the standard 
three-band copper-oxygen model of high- Tc cuprate su- 
perconductors. The model contains the physics of Zhang- 
Rice singlets at low doping and strong correlations, but 
for higher doping or weaker correlations the behavior be- 
comes band-like. We have examined the model both at 
strong coupling, where the state at half filling is a charge- 
transfer insulator-, and the physics at low doping is of 
Zhang-Rice singlets, and at weaker correlation strength, 
where antiferromagnetism is needed to make the model 
insulating at half filling and the appropriateness of the 
Zhang-Rice picture is less clear. We find that in both 
cases the model reproduces the low doping behvior but 
qualitatively fails to reproduce the experimental data at 
higher doping.— We conclude that if the experiment in- 
deed measures the oxygen density of states without other 
complicating effects, then the data indicate not just a 
failure of the Zhang-Rice model but a breakdown of the 
entire three-band model picture, with some additional 
orbital, not included in the three-band model, becoming 
important. 

An additional consequence of our work is that by com- 
paring the energies of the "Zhang-Rice band" and "up- 
per Hubbard band" features to our calculations we are 
able to estimate the correlation strength of the materials. 
The comparison places the materials on the metallic side 
of the metal/charge-transfer-insulater phase diagram, in 
agreement with previous work.^Siii 

The rest of this paper is organized as follows: Section 



im describes the three-band model to be studied and the 
methods. In section Hill we present the numerical results. 
Section HV] is a conclusion and discussion. 



II. MODEL AND METHODS 

The three-band Emery model considers Cu 3c?a;2_j,2 
and O 2pa orbitals) ^^'^'^ which hybridize via a copper- 
oxygen hopping tpd- We shall also consider oxygen- 
oxygen hoppings, adop ting the form implied by the con- 
siderations of Ref. The model exhibits both a low- 
doping antiferromagnetic phase and a high-doping para- 
magnetic phase but because the experimentally interest- 
ing behavior occurs at high dopings, we restrict our atten- 
tion in the paramagnetic phase. We distinguish the oxy- 
gen sites displaced from the Cu in the x and y directions, 
adopt the basis = {dka,Px,ka,Py,ka)i introduce the 
copper-oxygen hopping tpd = l.GeV and oxygen-oxygen 
hopping tpp in the form proposed by Andersenji^ we com- 
pare tpp = l.leV and tpp = The interaction part of 
the Hamiltonian is 

Hint ^ U^Ud^Hdl, (1) 
i 

and in the calculation presented here we take U — 9eV 
(the precise value of U is not important as long as it is 
larger than about 5eV). The band theoretic part of the 
Hamiltonian may be represented as a 3 x 3 matrix: 
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2itpd sin !f 



Hsband = -2itpd sin Ep -I- 2tpp{cos kx ~ I) 
-2itpd sin 4tppSin^sin^ 



2itpd sin -f 
4tpp sin ^ sin ^ 



(2) 



The model has two important parameters: U and 
A — Ep — Ed- The physically relevant case is large U, 
in which case the model is in the charge-transfer regime^ 
and the physics is controlled by the ratio of A to the 
copper-oxygen hopping tpd- At "half filling" (one hole 
in the d-p complex) the single-site dynamical mean-field 
approximation predicts that the model exhibits a param- 
agnetic insulating phase for small A whereas the param- 
agnetic phase is metallic at large A. In previous work^^ 
we have determined that for U — 9eV, A = 2eV is on 
the insulating side, but not far from the metal-insulator 
phase boundary, whereas A = AeV lies on the metal- 
lic side of the phase diagram, but also not far from 
the boundary. The correct value of A for the cuprates 
is controversialiS; we therefore consider A-values corre- 
sponding both to paramagnetic metal and to paramag- 
netic insulating phases at half filling. 



The model is solved using the single-site dynamical 
mean-field approximatioUj^'^ in which the key approxi- 
mation is a momentum-independent self-energy: 

/ E(z) \ 
E(z,k)^S(z)= 0, (3) 

\ y 

where z indicates real or Matsubara frequencies. The 
Green's function is 

G(Z, k) = {zl+fi~ S(z) - Hsband)"' , (4) 

and the spectral functions (density of states) of d and 
p orbitals are obtained from the imaginary part of the 
diagonal elements of the Green's function matrix. 

We employ two impurity solvers: Exact Diagonal- 
izations (ED)^^'^^ and the hybridization-expansion 
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continuous-time quantum Monte Carlo method 
(QMC)i2i. Because the two methods involve differ- 
ent approximations, comparison of the results helps 
confirm the accuracy of the methods. Analytic continu- 
ation of the QMC data is performed using the method 
described in Ref. The model has been previously 
studiedjiiii^Ti^i At strong correlations and low dopings 
the approximation reveals the "Zhang-Rice singlet" 
behavior with doped holes residing on the oxygen but 
strongly antiferromagnetically coupled to spins on the 
copper sites. As the doping is increased or the effective 
correlation strength decreased, the behavior reverts to 
moderately-correlated band behavior and a Zhang-Rice 
picture becomes inappropriate. 

From the calculation we obtain the electron spectral 
function A{u!) — ImG(aj) which we use to model the X- 
ray absorption experiment. As noted in Ref. !5:, the X- 
ray absorption cross section is not simply proportional 
to the product of a matrix element and electron spectral 
function because the hole in the oxygen Is core state 
interacts with the excited electron. We assume, following 
Ref. d that this effect provides a constant (Hartree) shift 
of the spectrum, and in particular has a negligible effect 
on the lineshape and on spectral weight. Therefore the 
x-ray absorption spectrum is given by 

B{uj) = C ■ Ao^y{uJ - cjo)[l - /(w - cjo)], (5) 

where Aoxy{^) is the calculated electron spectral func- 
tion, projected onto the oxygen site, [1 — f{uj)] is the 
complement of the Fermi function, restricting the result 
to unoccupied states, C encodes the matrix element and 
Wo is the energy difference between the final and initial 
state, including the effect of the excitonic interaction be- 
tween the core hole and the excited electron. In our cal- 
culations we set ojo for 0.1 hole doping case such that 
the upper Hubbard band peaks at 530.26^.4 The Wq for 
other doping values are computed by adding a Hartree 
shift to the 0.1 hole-doped case, using the Hamiltonian 
and parameters described by Ref. i5| with the calculated 
occupation numbers. 
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FIG. 2: Unoccupied portion of three-band model many-body 
density of states, projected onto oxygen orbitals, calculated 
as described in the text for dopings x indicated on figure us- 
ing both analytical continuation of QMC [panels (a)] and ED 
[panels (b)] impurity solvers for two values of charge trans- 
fer gap A; A = lieV (upper panels) corresponding to half 
filled paramagnetic metal and A = 2eV (lower panels) cor- 
responding to paramagnetic insulating phase at half filling. 
Parameters: U = 9eV , tpp — 0, tpd = 1.6eV and temperature 
T = O.leV. 



III. RESULTS 

Fig. [2] shows the doping dependence of the electron 
spectral functions, calculated as described above and pro- 
jected onto the oxygen orbitals, for A — 2eV (insulator 
in undoped case) and AeV (metal in undoped case), in 
absence of tpp. Results obtained from both analytic con- 
tinuation of QMC and ED impurity solvers are shown; 
the close correspondence between the results of the two 
methods indicates that the solution of the model is reli- 
able. 

Two features are evident in the calculated spectra: a 
higher energy feature which corresponds to the upper 
Hubbard band (its weight is small because this band is 
comprised mainly of d-states so it has a small projec- 
tion onto the oxygen states) and a lower energy feature 



corresponding to the Zhang- Rice band. The difference in 
energy between the two features is an estimate of the cor- 
relation gap in the system. One sees that the gap size is 
about AeV in the A = 2eV (insulator in undoped case) 
case, and around 2eV in A — AeV (metal in undoped 
case) case, and we expect the gap will further shrink if 
we go deeper in the metallic phase of the paramagnetic 
phase diagram. Experiments^'^ show that the observed 
gap between the two peaks is < 2eV , which is more con- 
sistent with the A = 4ey results, in agreement with our 
previously published papers^"' placing the cuprates on 
the paramagnetic metal side of the phase diagram, so 
that the correlations are intermediate rather than strong. 

Fig. [3] shows results for the case tpp — l.leV obtained 
from analytic continuation of the QMC data for the un- 
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FIG. 3: Calculated unoccupied oxygen density of states com- 
puted as a function of doping as described in the text for 
non- vanishing oxygen-oxygen hopping tpp = l.leV by ana- 
lytic continuation of QMC calculations. Other parameters: 
U = 9eV, tpd — l.GeV and temperature T — O.leV. 
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FIG. 4: Calculated relative intensity of the peak associated 
with the upper Hubbard band of O K-edge x-ray absorption 
spectra v.s. doping. The results are normalized such that 
their values at 0.22 hole doping are the same. The dashed 
line is only guide to eyes. Parameters: U = 9eV, tpd — l.GeV , 
T = O.leV, tpp and A = — £d as shown in the legends. 



occupied oxygen density of states. The similarity of these 
curves to the resuhs shown in Fig[5]indicates that oxygen- 
oxygen hopping does not have an important effect on the 
results. 

We have integrated the area in the lower energy feature 
shown in Figs. [2] and [31 normalized the results to the value 
at hole doping 0.22 and have plotted the resulting spec- 
tral weight as diamonds etc. in Fig.[T] For both values of 
A the integrated spectral weight is found to increase lin- 
early with doping over a wide doping range. Differences 
between the A = 2eV case (strong correlation, doped 



charge-transfer insulator) and the A = 4eV case (inter- 
mediate correlation, charge-transfer metal) appear only 
at low dopings x < 0.11. For completeness we present 
also in Fig. 2] the integrated oxygen weight in the upper 
Hubbard band. We see that the spectral weight of the 
Zhang-Rice band continue to increase at the overdoped 
region, while the weight of the upper Hubbard band de- 
creases. Both the weight in the Zhang- Rice band and the 
weight in the upper Hubbard band vary smoothly with 
doping for all parameters. The sharp break in the data is 
not observed in the calculation for any choice of parame- 
ters. We therefore conclude that the experimental paper 
understated the significance of the results: the data in- 
dicate not just a breakdown of the Zhang-Rice picture, 
but a failure of the three-band model itself. 



IV. CONCLUSION 

In this paper we have used single-site dynamical mean- 
field theory to study the Emery three-band copper-oxide 
model related to the cuprates. We have calculated the 
doping dependence of the intensity of oxygen K-edge x- 
ray absorption spectra. At high doping, our calculations 
does not reproduce the results of a recent experiment 
(Ref. 0). This implies a breakdown of "Zhang- Rice sin- 
glet" approximation and even a failure of the entire three- 
band model picture. 

We may speculate on the reason for the failure of the 
three-band model in the overdoped regime. One possibil- 
ity is that the electronic structure changes in such a way 
that new degrees of freedom (beyond the conventionally 
studied copper d^i^yi and oxygen po.) states become im- 
portant (for example other states in the d-multiplet or 
apical oxygen states) providing a new channel for adding 
doped holes. An alternative possibility might be that 
for some reason the arguments relating the experimen- 
tal measurement to the oxygen density of states break 
down, or that there is an unusual change in the matrix 
element. Such a change would presumably be related to 
a change in the electronic structure. A third possibility 
is that additional physics, such as a Upp on the oxygen 
states, begins to play a role; however it is not clear why 
this would lead to a sudden change in doping dependence. 
The effect of additional orbitals is presently under inves- 
tigation. 
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